A Theory of Ferroelectric Phase Transition in SrTiOs induced by 

Isotope Replacement 

Yasusada Yamada^-*, Norikazu Todoroki^^ and Seiji Miyashita^-* 
l)Advanced Research Institute for Science and Engineering, Waseda University 3-4-1, Okuho, 

Shinjuku-ku, Tokyo 169-0072, Japan 
2)Department of Applied Physics, The Tokyo University, Bunkyo-ku, Hongo 7-3-1, Tokyo 

113-8656 
(February 2, 2008) 

Abstract 

A theory to describe the dielectric anomaUes and the ferroelectric phase 
transition induced by oxygen isotope replacement in SrTiOa is developed. 
The proposed model gives consistent explanation between apparently contra- 
dictory experimental results on macroscopic dielectric measurements versus 
microscopic lattice dynamical measurements by neutron scattering studies. 
The essential feature is described by a 3-state quantum order-disorder sys- 
tem characterizing the degenerated excited states in addition to the ground 
state of TiOg cluster. The effect of isotope replacement is taken into account 
through the tunneling frequency between the excited states. The dielectric 
properties are analyzed by the mean field approximation (MFA), which gives 
qualitative agreements with experimental results throughout full range of the 
isotope concentration. The phase diagram in the temperature-tunneling fre- 
quency coordinate is studied by a QMC method to confirm the qualitative 
validity of the MFA analysis. 



I. INTRODUCTION 



Recently, the giant anomalous dielectricity and induced ferroelectricity by isotope re- 
placement of ^^O for ^^O in SrTiOa have attracted much attention from the view point of 
quantum ferroelectricity. [1-4] In 1999, Ito et al. [1] first reported the appearance of broad 
maximum of the dielectricity at T = 23K in SrTi^^Os suggesting the stabihzation of fer- 
roelectric phase upon isotope replacement of ^^0 to ^^O. Later, the experiments have been 
refined by using single domain specimens and also using partially doped SrTi(^^0i_a;^*0a;)3 
(hereafter ST018-x) materials. 

The results by Wang and Ito [3] on the dielectric constants are reproduced in Fig. 1 in 
somewhat modified form using a three dimensional representation. As is clearly observed in 
the figure, the peak of static dielectric constant at T = increases as the ^^O concentration 
X is increased to xc — 33%. Upon further increasing x, the peak position Tc{x) shifts 
continuously, and reaches Tc = 21 AK at x = 84%. When projected on x — T plane, the 
curve connecting the peak positions (the thick solid line in the figure) effectively represents 
the para-ferroelectric phase diagram. 

Ito et al. analyzed the experimental data of STOlS-x based on the 'multi component 
vector model' by Schneider et al., [5-7] which deals with the anharmonic lattice vibrational 
system with emphasis on the quantum effect manifested in zero-point vibration. In their 
analysis, the observed x — T phase diagram is reproduced if the 'critical quantum displacive 
hmit' (see Sec. 2) is assumed to be realized x — xc{— 33%). The temperature de- 
pendences of the dielectric constant e{T]x) with various values of x are analyzed based on 
Barrett formula [8,9]: 

in the temperature region of 30K < T < lOOK, taking Ti and Tq as adjustable parameters, 
and C is a constant. It is concluded that when Ti is varied in the range of 71K < Ti < 87K 
while To kept constant, the results show satisfactory agreements throughout the observed 
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x- values. Thus, as far as the macroscopic dielectric characteristics such as phase diagram 
and static dielectric constant are concerned, the experimental results seem to be explained 
satisfactorily by the anharmonic lattice vibrational system with quantum effect taken into 
account explicitly. 

Recently, however, microscopic observations of the soft TO phonon have been carried out 
on STOlS-x using neutron inelastic scattering technique. [10,11] The results show essential 
discrepancies with the expected characteristics of the TO phonon frequency u;o{T). The 
results are summarized in Fig. 2. In the figure the calculated curves of uJo{T) for an 
anharmonic phonon system (see Sec. 2) are plotted for the cases of a; = 0%, x = xc and 
X — 100%, together with the observed points Bit x — 0%, 33%, and 89%. As is clearly seen in 
the figure, the observed values for x — 33% and 89% are in complete disagreement with the 
calculated curves. That is, while the calculated curves give the remarkable characteristic 
that they reach the stability limit (a; = 0) at T > 0, the observed points do not seem to 
exhibit such tendencies. In fact, the soft TO mode frequency is practically x-independent. 
These results seem to indicate that there should be some other freedom of motion which is 
responsible for the anomalous dielectric behavior of ST018-X. It is the purpose of the present 
investigation to present a possible model to give a consistent explanation for the apparent 
discrepancy between the macroscopic versus the microscopic observations. We propose a 
model where we introduce a variable specified by 'bipolarons' due to the ^^O substitution. 
The model is a generahzed Blume-Capel (BC) model [12] with quantum fluctuation. In 
this model, the new degree of freedom causes the order-disorder phase transition and also 
zero temperature fluctuation of polarization. We analyse this model by the mean fleld 
approximation (MFA) and Quantum Monte Carlo (QMC) simulation, and compare the 
results with the observed data. 
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II. QUANTUM DISPLACIVE SYSTEM VS QUANTUM ORDER-DISORDER 

SYSTEM 



Ferroelectric phase transition is conventionally divided into two categories: displacive 
type and order-disorder type. Correspondingly, the quantum effect of polarization fluctua- 
tion is described somewhat differently. In this section we briefly review the quantum effect 
in both cases. 



A. Quantum displacive system 

The Hamiltonian to describe the system is given as follows [13]: 

^ = ^ E(^' + ^LQ') + ^ E a' + E v^JQiQ,, (2) 

i i i,j 

where Pi and Qi are respectively the momentum and amplitude of the local optical vibra- 
tional mode at the i-th site, a;ioc is the characteristic frequency, and Vij is the interaction 
between the induced polarization on different sites. Through the renormalization process 
of the anharmonic part, the system is represented by a quasi-harmonic system with the 
renormalized frequency cuk, 

n = lT.iPk+^lQl), (3) 

k 

cu',=ujI + X{Q'), (4) 

^o=^ioc-n? = 0). (5) 

Here, ujq stands for the 'bare' harmonic phonon frequency of the uniform [q = 0) TO mode. 
Notice that in order to undergo ferroelectric phase transition loq should be negative. That 
is, without the anharmonicity the system is intrinsically unstable against uniform TO mode. 
The renormalized frequency cul is temperature dependent through 

(g^) = ^coth^, (6) 
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which is extrapolated at higher temperatures as 

iQ') - ¥ (7) 

^loc 

giving the 'classical limit'. The renormahzed frequency uq is nothing but the so-called 'soft' 
mode frequency characterizing the displacive ferroelectricity 

Since the dielectric susceptibility (or the dielectric constant s) due to the soft TO mode 
is given by 

Xph^ -Zo, (8) 

the overall characteristic feature is conveniently overviewed by plotting Xph temperature 
(Fig. 3 (a)). As is evident in the figure, the effect of the quantum fluctuation shifts the Xph 
from the dashed curve representing the classical limit to the solid curve. 



B. Quantum order-disorder system 

In order-disorder type systems, the local potential for the order parameter has the multi- 
minimum structure. Hence, the variable associated with the polarization is considered to 
take several discrete values which are effectively expressed by a 'spin' variable. [14,15] 

In the simplest case of the double minimum potential, the effective Hamiltonian is ex- 
pressed in terms of the 2x2 Pauli spin operators as follows, 

n = hQY.^t-j:Jv^tcTh (9) 

i i,j 

where fl is the tunneling frequency of the hopping motion of the polarization, and Jij is the 
interaction between the local polarizations at the sites i and j. 
Within MFA the spin susceptibility is given by 

Xsp = (i_J(0);^sj Xsp, (10) 

where 
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(11) 



Xsp 



kBT' 



Here, Xsp stands for the susceptibility of an isolated tunneling particle (single particle sus- 
ceptibility), while the term in the bracket in Eq. (11) is the 'enhancement factor' due to 
cooperativity. Notice the classical hmit of xL is given by 



In order to overview the characteristic feature, we again plot Xsp vs temperature (Fig. 3(b)). 
We obtain a similar feature to the case of displacive type except that the quantum effect is 
due to the tunneling motion rather than the zero-point vibration of the particle. 



As has been already pointed out by several authors, [16-18] in spite of apparently quite 
different formalism of Xph and Xsp (Eqs. (8) and (10)), both can be reduced to the common 
expression called Barrett formula (Eq. (1)) except that the characteristic parameters, Tq, 
Ti, and C, are expressed in terms of different physical quantities as listed in Table I. 

This situation implies that as far as the macroscopic characteristics such as static dielec- 
tric constant are concerned, it is difficult to distinguish which category the system under 
investigation belongs to. In order to distinguish the category, some microscopic measure- 
ments are necessary. In the case of SrTiOs, neutron inelastic scattering gave the decisive 
information. In 1968, Yamada and Shirane [18] observed the phonon dispersions of ST016 
by neutron inelastic scattering technique. They found the well-defined TO phonon branch 
which shows strong softening particularly at g = 0. The observed temperature variation of 
the soft mode frequency u!o{T) was fitted satisfactorily to Barrett formula in a wide tem- 
perature range of 20K <T < 300K. The results clearly indicate that ST016 belongs to the 
typical 'quantum displacive' system. 




1 



(12) 



kBT' 



C. Barrett formula 
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III. 3-STATE QUANTUM ORDER-DISORDER SYSTEM 



As a natural extension of ST016 described in the last section, one expects that ST018-X 
would be understood within the same framework of quantum displacive system. However, 
recent neutron scattering studies on ST018-33 [10] and ST018-89 [11] seem to provide es- 
sentially contradictory results to this anticipation. As is summarized in Fig. 2, the observed 
temperature dependences of the soft TO phonon frequencies are practically x-independent. 
On the other hand, the theoretical values based on the quantum displacive system. 



with To and Ti given by fitting to the observed e{T; x), show substantially different temper- 
ature dependences for the typical cases oi x — xc{— 33%) and x — 100%. 

It is instructive to plot at the zero temperature limit (T ^ 0) as a function of x as 
given in Fig. 4. It is seen that the contribution due to the zero-point motion of soft 
TO phonon which is deduced from the neutron data (the sohd squares and the dashed 
line in the figure) is irrespective of the observed critical behavior of £{T 0;x) around 
X = 33%. This fact suggests that the dielectric behavior of ST018-a; is composed of two 
different components, and that the anomalous part in low temperature region originates 
from the quantum fluctuations associated with some other freedom of motion than the soft 
TO phonon mode. 

Recently, Hasegawa et al. [19] reported the observation of a giant dielectricity induced 
by laser irradiation of ST016 crystals, which suggests the possibility that the anomalous 
dielectricity in SrTiOs is due to the polarization fluctuation associated with the electronic 
excited states. The excited state is created by single electron transfer from one of the oxygens 
to titanium: Ti^+Og^^Ti^+Og^O". Since the excited state breaks inversion symmetry 
inducing spontaneous distortion of the octahedron, the adiabatic potential for the polar 
local mode, Qioc, would have the characteristic as shown schematically in Fig. 5. It is 
noticed that in addition to the minimum at Qioc = 0, corresponding to the ground state. 




(13) 
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there are local minima at Qioc = ±Qo, corresponding to the excited states. The excited state 
accompanied by the spontaneous distortion may be viewed as a pair of an electron polaron 
around Ti^"*" and a hole polaron around 0~^. This pair is a kind of 'bipolaronic state'. 

Based on these considerations we construct a model which takes into account the bipo- 
laronic exited states in addition to the ground state of TiOe cluster. If we neglect the small 
amplitude vibrations around the local minima in Fig. 5, the freedom of motion of the system 
is simply represented by a 'spin' variable which takes three discrete values; — 0, +1, and 
—1. We further assume that while the transition between the ground state and the excited 
states 5^ = ±1 is allowed only by thermal excitation, quantum hopping (or tunneling) 
is allowed between the degenerated bipolaronic states {S^ — 1 <-> —1). 

The 'spin' part of the freedom of motion is characterized by a 3-state model as the BC 
model. In orderer to take into account the quantum hopping, we introduce a quantum 
transfer term. Thus, the Hamiltonian of the model is given by 



i i,j i 



z\2 



(14) 



Here, 5*? and are the spin operators expressed by 3 x 3 matrices rather than the 2x2 
Pauli spin operators a in Eq. (9): 



'^0 1^ 



Sf 





1 



^o 0^ 



Sf 




0-1 



(15) 



The interaction Jij denotes the ferroelectric interaction between the polarization. Here we 
consider interactions only between the nearest neighbour spins. It should be noted that 
is not the standard operator for the x-component of the angular momentum. The parameter 
A in the third term describes the excitation energy to create the bipolarons. Within MFA, 



the Hamiltonian is given in matrix form by, 

/ 



'Hmfk — 



-Jz{S') + A hVt 


hVL Jz{S') + A 



(16) 



where z denotes the numbers of the nearest neighbour sites. 

Through a standard process, the free energy of the system is exphcitly given by 



F(0 = -TlnTr (e~^^^^/^) + \-Jzi^ 



2 ^ 



Tin 



2 cosh ^^i^2 + ^2j2^2 + exp ^ 



T 



(17) 



where ^ = {S^) stands for the order parameter or the spontaneous polarization. For concise- 
ness of the expression we have put kB — h — 1. The quahtative nature of ordering studied 
by the MFA are confirmed by QMC simulation as given in the Appendix. 



IV. ANALYSIS OF EXPERIMENTAL RESULTS 

Once the general expression of F{^) is given, we can obtain various thermodynamical 
quantities such as dielectric constant e{T; x) as a function of the three parameters, ft, J, and 
A (essentially two independent parameters, Q/J and A/ J). In order to explain the experi- 
mental results, we have to determine the dependence of the parameters on the concentration 
X. We conjecture that among these parameters, the tunneling frequency Q would be most 
sensitive to the isotope replacement. Hereafter, we discuss only the limited case that only 
fl is x-dependent and A/ J — 1. The second assumption implies that the excitation energy 
to create a bipolaron is comparable to the total nearest neighbour interaction. 

It is easily seen that, in the model (16), the quantum mechanical stability limit for the 
ferroelectric state is 

n{xc)/J = 6.0 (18) 

for the case of 2; = 6 corresponding to the simple cubic lattice. We recognize that this 
condition is satisfied at a; = 33% (= xc ) in STOlS-x. A linear dependence of fl{x) 

n{x)/J =n{xc)/J + 0.9{xc-x), 0<x<l, (19) 

is adopted as the simplest choice. The parameter value 0.9 has been chosen so that overall 
phase diagram becomes qualitatively consistent with the observed one (see Fig.6(a)). This 
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dependence means that the tunnehng frequency changes by 15% within the full range of 
isotope replacement: 0% < x < 100%. 

The results of the calculations are summarized in Fig. 6. In Fig. 6(a), the ordered phase 
is specified as 'ferroelectric' within the quotation mark, which implies that the ordered state 
is somewhat different from the ordinary ferroelectric state. In the case of the ordinary 
ferroelectric states, the polarization fluctuation is substantially replaced by the spontaneous 
polarization as T — > 0. In contrast, in the present case there remains considerable amount 
of polarization fluctuation even at T = which is associated with the tunneling fluctuation 
of the spins as well as the soft mode vibration. In Fig. 6 (b) and (c), the 'experimental' 
values are given after subtracting the contribution due to the quantum displacive part given 
by Eq. (8) in order to facilitate direct comparison with the calculations based on the 3- 
state quantum order disorder model. Although the calculated results are still far from 
the situation to carry out quantitative analysis to fit the experimental data, the overall 
qualitative features seem to be satisfactorily reproduced. 

V. CONCLUSIONS AND DISCUSSIONS 

A theory to describe the dielectric anomalies induced by oxygen isotope replacement in 
SrTiOs is proposed which gives consistent interpretation between the macroscopic properties 
(static dielectric constant) and the microscopic properties (TO phonon frequency observed 
by inelastic neutron scattering) . The essential feature of the anomaly is characterized by the 
3-state quantum order-disorder system (14) associated with the bipolaronic excited states 
in addition to the ground state of TiOe cluster. By taking into account the effect of the 
isotope replacement through the frequency of the tunneling between the bipolaronic states, 
the dielectric properties are investigated. The results give qualitative agreements with the 
experimental results throughout full range of the isotope concentration, < x < 100%. 

In the present treatment, we have assumed for simplicity that there exists a single (non- 
degenerate) polar mode, Qioc, which couples to the electronic excitation. In the case of TiOe, 
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however, the polar modes are triply degenerated forming the basis functions, Qx-, Qy, Qz, oi 
Ti„ irreducible representation of the cubic point group m3m. Therefore, the adiabatic po- 
tential of the polar mode should have quasi-spherical symmetry in the {Qx, Qy, Qz) space. 
There are at least six minima located on the sphere with a fixed radius Qq and the lowest 
potential barriers (saddle points) between the minima through which the bipolaron under- 
goes tunneling are also situated on the same sphere. That is, the 'spin' is rotating along 
the spherical ditch as it changes the state (see Fig. 7). Correspondingly, the tunneling 
modes would be split into six levels instead of two as assumed in the present paper. In 
order to carry out quantitative analysis proper symmetry considerations would be crucially 
important. 

The direct evidence of the vahdity of the present picture should be given by the obser- 
vation of the tunneling excitations using Raman scattering, neutron scattering, etc.. Up 
to this date, any well-defined excitation spectrum has not been reported in addition to the 
ordinary phonon excitations observed by neutron scattering technique. Recently, Noda et 
al. [20] observed quasielastic neutron scatterings around the Bragg reflections which exhibit 
strong temperature dependence. Therefore it would be possible that the tunneling modes 
are heavily damped so that the scattering spectra due to tunneling fluctuation have resulted 
in so-called 'central peak' spectrum centered at a; = 0. Similar kind of critical over-damping 
of the soft mode is commonly observed at various structural phase transitions. More detailed 
spectroscopic observations of the critical fluctuations around x < xq Bi lower temperatures 
are needed. 

APPENDIX A: QUANTUM MONTE CARLO SIMULATION 

We study the phase diagram of the model (14) by a QMC method with the Suzuki- 
Trotter decomposition. In the present model non-diagonal interaction exists only at each 
site. Thus, the model is simply transformed to a four-dimensional classical model. We 
performed Monte Carlo simulation on this model by using a cluster heat bath algorithm. 
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[21] Recently efficiently of QMC method has been extensively developed by introducing the 
loop algorithm to avoid the slowing down due world line algorithm at low temperatures. 
[22] The loop algorithm has been also extended to system s with large S. [23] In this study, 
however, we adopt a cluster heat bath algorithm which provides an efficient sampling in 
strongly spacially anisotropic systems. In the cluster heat algorithm, the Boltzmann weight 
is calculated straightforwardly by using one dimensional transfer matrix method. This 
method is efficient for the quasi one-dimensional discrete spin models. In this method we 
do not need to provide the graph representation of the cluster which depends on the details 
of the model. Thus, this cluster heat algorithm algorithm is very convenient for the present 
study. 

The Suzuki- Trotter decomposition transforms the Hamiltonian (14) into 

N 



where TieS is a four-dimensional classical model with a variable a — —1,0, and 1 



(Al) 



i^iv = ^log(cot^), (A3) 



A 1, 



1 , f2(3n\ 

- smh 



(A5) 



2 \ N J 

where N is the Trotter number and we adopt the periodic boundary condition for the Trotter 
direction a^+i = (y\- 

We consider the /^-th chain along the Trotter direction whose Hamilton is given by 

JV JV _ JV 

i=\ 1=1 i=l 

where = JJ2(jj,i/) <^i,ij.- We define the following transfer matrix: 
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^ exp{X^((77v(7i + (7i(72 + (72(73 H h (7n_i(7 



+ {Hiai + i/2CX2 + • • • + i/n(7) 



(A7) 



Tn{aN, cr) can be readily obtained from the recursion formula 



7;((Tjv, (x) = X! cr) exp{Kan-ia + if„(T - A(t^}, (n > 2) (A8) 



with 7j((T7v,cr) = exp(_ft'Ar(TAr(T + -ffi(Ti — A(T^). The probability PNicr^ = 0, ±1) that the 



boundary spin takes 0, ±1 is given as 




(A9) 



We can determine the value of ctat by using this probability. The other spins ak{k — N — 
1, N — 2, ■ ■ ■ ,1) are determined by the probability given by 



Thus, we determine new spin configuration of the chain. 

In Fig. 8(a) we present the phase diagram in the coordinate (Q, T) obtained by the QMC 
simulation. In Fig. 8(b) the phase diagram obtained by MFA is also given, We confirmed 
that the qualitative feature of the phase diagrams is common in these two cases. Here it is 
to notify that a ferromagnetic phase appears by being induced by quantum fluctuation in 
the case of large A region (A ~ 3). That is, at A = 3 the ferromagnetic phase is suppressed 
in classical BC model (0 = 0), but it appears at finite value of fl. 



(AlO) 



E(T^v((7Ar, (7) ' 



where 



Tk{aN,cr) = 7;exp(X(7(7fe+i). 



(All) 
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FIGURES 




FIG. 1. (Color) Temperature dependences of the static dielectric constant of ST018-X at vari- 
ous x-values observed by Wang and Ito. [3] The thick solid line on the T-x plane (top of the graph) 
indicates the phase boundary between the paraelectric phase and the ferroelectric phase. 
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FIG. 2. Comparison of the experimental soft TO-mode frequencies ujq(T) observed by neutron 
inelastic scattering and the theoretical curves based on the anharmonic quantum lattice vibrational 
system for the cases of x = 0% (dashed carve), at the quantum displacive limit x = xc (solid curve), 
and X = 100% (dotted carve). 



19 




(a) (b) 
FIG. 3. Schematic representations of the effect of quantum fluctuations appearing in ~ T 

curves, (a) The susceptibitity Xph (Eq. 8) for the case of quantum displacive system where the 

quantum effect is due to the zero-point vibration, and (b) Xrmsp (Eq. 10) for the case of quantum 

order-disorder system where quantum effect is due to the tunnehng of the quasi-spin variable. The 

dashed lines correspond to the cases for classical limits. The effects of quantum fluctuations are 

visualized as the differences between the dashed lines and the solid curves. 
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X 

FIG. 4. The observed dielectric constants as extrapolated to T = 0. The solid circles are the 
values by macroscopic dielectric measurements. The solid squares are the values deduced from the 
microscopic measurements of ujq{T) using (8). The solid and dashed lines are for guides of eyes. 
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-go Qo " Sloe 

FIG. 5. Schematic diagram of the adiabatic deformation potential of TiOg cluster as the ampli- 
tude of the local polar mode is varied. The central valley is associated with the electronic ground 
state, while the side valleys are associated with the first excited state, which is characterized as 
the bipolaronic state. 
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FIG. 6. Comparison of the observed characteristic features in e{T;x) (left column) and the 
corresponding theoretical curves (right column) with parameter values given in the text, (a)phase 
diagram QP: quantum paraelectric phase 'F': extraordinary ferroelectric phase (See the text.) 
(b) Temperature and the concentration dependence of £{T;x) In the experimental curves the con- 
tribution of the TO soft mode fluctuation has been subtracted. (c)Dependence of £{T 0;x) 
on the concentration. The experimental values are given by subtracting the contribution due to 
quantum displacive fluctuations. 
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FIG. 7. The adiabatic potential for the triply degenerated polar modes {Qx,Qy,Qz) within 
(Qx, Qy, 0)-plane. The minima are at (ibQojO, 0) and (0, ibQo,0), and the lowest barriers are at 
(±, Qo, ±Qo) 0). The tunneling bipolaron is 'rotating' along the ditch as it changes the state. 
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(a) 



(b) 



FIG. 8. Phase diagram in the cordinate (O, T) by using QMC simulation (a) and MFA (b) 
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